Formation of Nanopillar Arrays in Ultrathin Viscous Films: 
The Critical Role of Thermocapillary Stresses 



o 
o 

< 



c : 
i : 

q: 

O 1 

• i-H ' 

C/5 ■ 
>Y 

^ : 

Ph. 



> 

&\ 
OO 

o 

OS 

o 



Mathias Dietzel and Sandra M. Troian * 
California Institute of Technology, 
1200 E. California Blvd. MC 128-95, Pasadena, CA 91125 

Experiments by several groups during the past decade have shown that a molten polymer nanofilm 
subject to a large transverse thermal gradient undergoes spontaneous formation of periodic nanopil- 
lar arrays. The prevailing explanation is that coherent reflections of acoustic phonons within the film 
cause a periodic modulation of the radiation pressure which enhances pillar growth. By exploring a 
deformational instability of particular relevance to nanofilms, we demonstrate that thermocapillary 
forces play a crucial role in the formation process. Analytic and numerical predictions show good 
agreement with the pillar spacings obtained in experiment. Simulations of the interface equation 
further determine the rate of pillar growth of importance to technological applications. 

PACS numbers: 47.85.md,63.22.-m,47.20.Dr,47.61.-k 



Fabrication of high resolution, large area arrays for 
micro-optic, photonic and optoelectronic devices relies 
heavily on photolithographic patterning by projection of 
UV light. This process is inherently slow and costly due 
to multiple step-and-repeat procedures required for de- 
position, exposure and removal of photoresist layers for 
construction of 3D components. Harsh developer and 
etching solutions also imbue structures with significant 
surface roughness, which ultimately limits performance 
due to scattering losses. Photoresist masks must also be 
UV compatible, which restricts the assortment of ma- 
terials used in processing. Alternative patterning tech- 
niques based on "resistless lithography" may usher a new 
era in lithographic patterning in which structure forma- 
tion is obtained by directed deposition or mechanical em- 
bossing of material as with microcontact printing, mi- 
cromolding, microembossing and nanoimprinting [l], 0] • 
An even newer approach for emergent technologies re- 
lies on film patterning by hydrodynamic instabilities in 
nanoscale films as with templated dewetting 0] . The use 
of fluid instabilities for controlled formation of large area 
patterning provides an interesting route for non-contact 
fabrication of periodic microscale and nanoscale arrays. 

Several groups pursuing this approach have inves- 
tigated the spontaneous formation of pillar arrays in 
molten polymer nanofilms subject to a large transverse 
thermalgradient [1, [1,0 SB Si- The prevailing explana- 
tion [7|, [8| is that coherent reflections of acoustic phonons 
(AP) within the film cause periodic modulation of the 
radiation pressure, similar to an acoustic Casimir inter- 
action force This normal stress opposes capillary 
forces and enhances the growth of protrusions. Such a 
process, however, requires coherent phonon propagation 
within a molten amorphous polymer film and an aver- 
age phonon mean free path at least as large as the film 



thickness. Experiments have shown that solid nanoscale 
polymer films at temperatures — 193°C < T < 27°C can 
support phonon attenuation lengths of about O(IC) 1 — 10 2 ) 
nm at frequencies in the 100 GHz range [ll[. However, 
such long propagation lengths have never been demon- 
strated and are considered unlikely in amorphous molten 
films (glass transition temperature 100 - 120 °C). By 
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FIG. 1: (a) Experimental setup for formation of nanopillar arrays, 
(b) AFM image of PMMA pillars [j|: d = 260 nm, h a = 95 nm, 
Amax = 3.4^im, AT unknown, (c) Optical micrograph of PS pillars 
[H: d = 345 nm, h a = 100 nm, A max = 4.5/mi, AT = 46 °C. 
(d) AFM image of PMMA pillars 0: d = 163 nm, h a = 100 nm, 
A max = 6.5Mm, AT = 10 °C . 

investigating an unexplored limit of interfacial instability 
easily accessible to nanoscale films, we demonstrate that 
thermocapillary (TC) forces play a crucial if not dom- 
inant role in this formation process. According to this 
mechanism, perturbations in film thickness generate pe- 
riodic disturbances in the surface temperature which lead 
to periodic modulation of the interfacial thermocapillary 
stress. These tangential stresses cause the polymer melt 
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to flow toward cooler regions, thereby driving protrusions 
toward a cooler target substrate. These elongated pillars 
are observed to solidify rapidly in place as soon as the 
driving force is removed. Nanostructures fabricated in 
such non-contact fashion from the molten state are ex- 
pected to exhibit superior optical performance due to the 
specularly smooth interface obtained upon solidification. 
Investigation of this new regime of instability is there- 
fore pertinent to fundamental studies of thcrmocapillary 
flow as well as technological innovations for non-contact 
resistless lithography. 

A schematic diagram of the typical experimental setup 
is shown in Fig. [Ha). Polymers like polystyrene (PS) 
or poly(methyl methacrylate) (PMMA) are dissolved in 
solvent and spun cast onto a silicon wafer to a thickness 
80 nm < h a < 130 nm. The coated wafer is overlaid by a 
second wafer and vertical spacers are used to maintain a 
separation distance d > h Q where 100 nm < d < 600 nm. 
Ratios of d/h Q in experiment range from approximately 
2-8. The bottom wafer is placed on a hot substrate held 
at temperature 130°C < T 2 < 170°C; the polymer free 
surface is cooled from above by proximity to the cold 
substrate held at T\ < T%. Both temperatures are main- 
tained above the polymer glass transition to ensure a 
molten film during the formation process. Despite the 
fact that 10°C < T 2 - Ti < 55°C is not large, the 
small gap size establishes a very large transverse gra- 
dient AT/d ~ 10 6 - 10 8 °C/cm. In the experiments 
under study, films were subjected to a thermal gradient 
overnight and then quenched to room temperature to so- 
lidify the structures formed. The top wafer was then 
removed and optical microscopy or atomic force micro- 
scope (AFM) images obtained, which revealed the pat- 
terns shown in Fig. [1] (b)-(d). 

It is well known that much thicker liquid films (cm to 
mm) subject to much smaller thermal gradients can de- 
velop periodic cellular patterns through Rayleigh-Benard 
(RB) or Benard-Marangoni (BM) instability [jj]. These 
instabilities, however, generate shallow corrugations not 
needle-like protrusions. Onset of instability requires that 
the critical Rayleigh number Ra c for buoyancy driven 
flow (which scales as h*) or the critical Marangoni num- 
ber Ma c for thcrmocapillary flow (which scales as h 2 ,) 
exceeds 660 - 1700 or 50-80, respectively, depending on 
boundary conditions. For the nanoscale films shown 
in Fig. [TI Ra w 10" 16 and Ma w 10~ 8 , ruling out 
these mechanisms as possible causes for pillar forma- 
tion. A lesser known deformational instability leading to 
more pronounced dry spots or elevations [l3| has recently 
been observed in microscale films (50 < h a < 250/z.m) 
14, Til, IH in which the thermocapillary stress is coun- 
terbalanced by capillary and gravitational forces. On- 
set of instability requires that the inverse dynamic Bond 
number D c = 7<rATfii m / pgh 2 > 2/3(1 + F)~ l , where p 
is the liquid density, jt = \dy/dT\, 7 is the liquid sur- 
face tension, ATfri m is the temperature drop across the 



liquid layer, F = (1 — k)/{D + k — 1) is an order one con- 
stant, D = d/h , and k = fc a h /fciiq is the ratio of thermal 
conductivities. Parameter values for the experiments in 
Fig. Q] reveal that D c > O(10 7 ) and G ~ O(10~ 14 ), far 
beyond regimes of instability previously investigated in 
which D c ~ O(10 _1 - 1) and G ~ O(10 _1 - 10 2 ). 



In this paper, we investigate unexplored consequences 
of this deformational instability for nanoscale films sub- 
ject to a very large transverse gradient. The govern- 
ing geometric and dynamic parameter ranges are con- 
strained by e 2 = (ZioMmax) 2 < 1, e.Re 0, eRePr -> 
and G — > i.e. gravitational stabilization is absent. 
The lateral scale for the slender gap ratio, e, is set by 
^maxj the wavelength of the maximally unstable mode 
(i.e. fastest growing mode) obtained from linear stability 
analysis, which in experiment corresponds to the average 
pillar spacing. Here, Pr — v/a is the Prandtl number, 
G = gho^/va is the Galileo number, Re = u c h /u is 
the Reynolds number, u c is the characteristic lateral flow 
speed set by thermocapillary forces, g is the gravitational 
acceleration constant, v is the polymer kinematic viscos- 
ity and a is the polymer thermal diffusivity, both evalu- 
ated at the temperature T 2 . These constraints establish 
the slender gap approximation for momentum and ther- 
mal transfer [171 ] in the limit where there is no hydrostatic 
restoring force. 



Within this approximation, the differential equations 
for momentum and energy conservation decouple com- 
pletely [l8j]. The energy equation reduces to a ID 
steady state process for thermal conduction across an 
air/polymer bilayer with an internal undulatory inter- 
face H(X,Y,t); variation of the interfacial temperature 
with time arises solely through displacement of this in- 
terface. The relevant variables are normalized according 
to (X,Y) = (x/Xl^y/Xl^), Z = z/h a , H = h/h , 
9 = (T - Ti)/AT where AT = T 2 - T t , t = u c t/L, 
P = eh p/r)u c and T = jr\u c where t, p and rj = 77^2) 
denote real time, pressure and melt viscosity. The conti- 
nuity equation for incompressible flow yields the relevant 
scaling for the velocity fields; namely U — (U, V, W) — 
{u/u c , v/u c , w/eu c ). The boundary conditions for the ve- 
locity and stress fields are the usual no-slip condition, 
impenetrability at the solid surface Z — 0, and a jump in 
the normal and tangential stresses at Z = H(X, Y, r). In 
the slender gap approximation, these stress conditions re- 
duce to P z =h = -C^ _1 V 2 i7 and V S T = -M~a\ s <d z =H 
where V s Q z =h = -kDV\\H/[D+{k-1)H]- 2 , the ther- 
mal Marangoni number Ma — ejT AT '/ '{rju c ) , the sur- 
face gradient operator V s — * V|| = (d/dX,d/dY) [lij ]. 
and Ca — rju c / (-ye 3 ). The kinematic boundary condi- 
tion, dH/dr = Wz=h, is then re-expressed to yield the 
4th-order, non-linear equation for the evolution of the 



3 



air/polymer interface 17] : 
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In our simulations, u c was set by the choice that the 
interfacial thermocapillary stress and V||£f equal unity 
i.e. {8U/dZ)z= H = dT/dX = L_such that u c = 
(47r) 2 7e 3 /3r? 0. With this choice, Ca = (4tt) 2 /3 and 
Ma" = 37 T AT/[(47r) 2 7 e 2 ] = [£) + (« — 1)] 2 /(kT>) - 

o(io- 1 - 10 1 ). 

For the parameters values pertinent to experiment, 
spatial gradients in iJ for times r < 1 are 0(1) ; conse- 
quently, the first term in Eq. ([1]) is 0(1), the second term 
is O(10°-10 1 ), and the third term is O(10" 2 ). The desta- 
bilizing thermocapillary term is therefore O(10 2 — 10 3 ) 
larger than the stabilizing capillary term, which is the 
reason that pillar-like structures can form in this system. 
A similar derivation leading to Eq. ([TJ) yields the corre- 
sponding equation for the acoustic phonon model (tJ: 
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H 3 \ 
-Vn-ff + -=Vf l H = 0, 
[k-1)H] 2 11 3Ca 11 J 

where Q = 2Qk a AT/(u p 'ye 2 ), Q is a phenomenological 
reflectivity coefficient, u p denotes the speed of sound in 
polymer, and e = ^o/A^ x . The liquid flux in Eq. ([I]) 
due to shear flow by thermocapillary forces is instead 
replaced in Eq. ([2|) by a pressure driven flow due to 
acoustic phonons reflecting from the film interfaces. 

Conventional linear stability analysis yields the (di- 
mensional) wavelength for the maximally unstable mode 
corresponding to these different models, namely 
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These characteristic length scales can be directly com- 
pared to the pillar spacings observed in experiment. Lin- 
ear stability also yields the dimensionless cutoff wave 
number K c below which all modes are unstable (dis- 
cussed further in the summary) . For the TC model, 
Kj c = (3/2MMCaKD)^ 2 /[D + (k - 1)]. Solutions of 
Eq. ([3]) shown in Fig. [2^a) contain no adjustable param- 
eters. While the functional dependence of Eq. ([3]) on d is 
in good agreement with experiment, the TC model sys- 
tematically over predicts the average pillar spacing, in 
some cases by as much as 35 - 40%. To further test the 
dependence of A™ x on d, we performed a least-squares fit 
of the experimental data to the form of Eq. |3j), namely 
dd 1 ' 2 + C 2 d- 1 / 2 , as shown in Fig.^b). 

Superimposed on these curves is also a least-squares fit 
to Eq. Q with fitting parameters Q = 6.2 and u p = 1850 
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FIG. 2: (a) Experimental data from Refs. [10, SHI- Material 
constants, evaluated at the temperature T2 = 170 °C, were ob- 
tained from Refs. [27], HE EE H3l ■ ( b ) Fitting coefficients are of the 
form [Expt A-D, C\ (10 3 M"i) ' 5 , C 2 (0.1/im) 1 - 5 ]: [A, 0.353, -34.7], 
[B, 0.650, -64.6], [C, 0.379, -46.0], [D, 0.340, -30.7]. For the AP 
model, Q = 6.2 and u p = 1850 m/s. 
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FIG. 3: Solutions of Eq.Ufor = 11 °C and 46 °C. 

m/s 0- We have determined that several experimental 
factors may contribute to the offset observed in Fig.[2ja); 
here we focus on a key issue involving the measured val- 
ues of h in the literature. Spin cast polymer films are 
prone to significant solvent retention 2lJ, [22J , typically 
eliminated by post annealing for several hours in a vac- 
uum oven at elevated temperatures. Significant shrinkage 
in film thickness has been reported based on ambient va- 
por pressure, time and temperature of the bake. In the 
experiments of Ref. @, d, @, , there are no reports of 
film annealing following spin casting, which would lead 
to overestimates of h . Fig.[3]indicates the strong depen- 
dence of A™ x on h and d and a sharp drop in the pre- 
dicted wavelength for small values of AT and h - While 
smaller values of h due to film shrinkage leads to very 
good agreement, we note that solvent evaporation is as- 
sumed to occur prior to insertion in the experimental 
assembly and therefore plays no role in pillar formation. 

To better understand this instability, we examined the 
Lyapunov free energy [24], [25| F — J £ dXdY for the 
evolving film where 



£=(V||#) 



3nMa Ca 
D 



H 



1 + X H 



+ln(l + x) 



(5) 

and x = {k — l)/D. The 1st and 2nd terms on the 
right represent the capillary and thermocapillary contri- 
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butions, respectively. Finite element simulations with 
periodic boundary conditions and 2nd order Lagrangian 
shape functions for spatial discretization of the film 
height were conducted. As shown in Fig. Ufa), once r 
exceeds unity, film evolution enters the nonlinear regime, 
as indicated by the steep drop in free energy due to the 
overwhelming influence of thcrmocapillary forces. The 
terminating points correspond to contact with the upper 
plate. This instability is non-saturating (i.e. no steady 
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FIG. 4: (a) Evolution of the Lyapunov free energy for the case 
h o =100 nm, d=285 nm, AT = 46°C. Real time conversion t = 
2.28r hrs. Inset: Plot of H(X,Y,t = 1.1025) from finite element 
simulations of Eq. JT). (b) Time required, ttop, for fastest growing 
pillars to contact upper wafer for parameter values given in (a). For 
variation with h , d=285 nm; for variation with d, /i o =100 nm. 

state solution in contrast to RB or BM instabilities); 
pillars will simply continue to grow until contact with 
the upper plate. Shown in Fig. 0JU) are predictions of 
the time required for the fastest growing pillars to con- 
tact the cooler plate. The TC curve was obtained from 
Eq.([T]); the AP curve is an estimate based on the growth 
rate corresponding to Eq. ([4]) from linear stability anal- 
ysis. Our estimates in Fig. 0Jb) a ls° indicate that the 
fastest growing nanopillars contact the upper substrate 
within a few hours of formation, while the experiments 
typically lasted overnight. Filaments bridging both sub- 
strates might then undergo capillary or even themocap- 
illary migration toward narrower gaps within the slight 
wedge geometry used in experiment (tilt rj 1/im per cm 
according to [23|). This secondary effect would also lead 
to better agreement due to smaller measured values of 
A™x. Much larger values of D = d/h a were also inves- 
tigated numerically. Evolution into the nonlinear regime 
corresponding to much larger pillar amplitudes generates 
in-plane hexagonal ordering, in agreement with experi- 
mental observations [1^1 . 

In summary, we have shown that thermocapillary 
stresses play a crucial if not dominant role in the for- 
mation of pillar arrays in molten nanofilms subject to a 
large transverse thermal gradient. The parameter range 
investigated corresponds to an unexplored limit of defor- 
mational instability in which destabilizing thermocapil- 
lary forces far outweigh stabilization by capillary or grav- 
itational forces. The predominance of thermocapillary 



effects allows the formation of elongated nanostructures. 
Our analysis indicates that any Newtonian liquid subject 
to even smaller thermal gradients will undergo pillar for- 
mation for wavenumber disturbances K < K c . However, 
technological demand for high resolution optical or other 
large area arrays with very small feature sizes suggests 
use of very large thermal gradients, smaller gap widths, 
and smaller film thicknesses h a . 
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